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We investigate the dynamics of two-dimensional site-diluted Ising antiferromagnets. In an exter- 
nal magnetic field these highly disordered magnetic systems have a domain structure which consists 
of fractal domains with sizes on a broad range of length scales. We focus on the dynamics of these 
systems during the relaxation from a long-range ordered initial state to the disordered fractal-domain 
state after applying an external magnetic field. The equilibrium state with applied field consists 
of fractal domains with a size distribution which follows a power law with an exponential cut-off. 
The dynamics of the system can be understood as a growth process of this fractal-domain state in 
such a way that the equilibrium distribution of domains develops during time. Following these ideas 
quantitatively we derive a simple description of the time dependence of the order parameter. The 
agreement with simulations is excellent. 
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> . I. INTRODUCTION 



One starting point for the theoretical investigation of the behavior of random magnets and related systems is to 
consider the system as consisting of finite ordered regions, called clusters, domains or droplets. In this article we 
want to consider the case of highly disordered systems where the magnetic structure consists of fractal domains with 
\Q ' domain sizes on a broad range of length scales and we will derive a simple theoretical description based on scaling 
ON ■ ideas. Nearly all of the assumptions we make and scaling relations we use will be proven in detail by exact ground 
' state calculations as well as by Monte Carlo simulations. 

The system we consider is the 2D diluted Ising antiferromagnet in an external magnetic field (DAFF) which can be 
used to study typical behavior of strongly disordered systems as there are metastability, slow dynamics, and domain 
structures M. Related systems like random- field systems which are thought to be in the same universality class like 
the DAFF jg] or spin glasses may behave similar. 
q The Hamiltonian of the DAFF in units of the next-neighbor coupling constant J reads 
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H = ( ,< ,n,n , - 2J BeiiJi (1) 
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C3 ' where Cj = ±1 denote Ising-spins and e,- = 0, 1 a quenched disorder. 

The phase diagram of the 2D DAFF consists of an antiferromagnetic low temperature phase for magnetic field B = 
and a disordered phase for all finite values of B. In contrast, in three dimensions there exists a long-range ordered 
phase also for small magnetic fields ||. For higher fields and low temperatures the DAFF develops a disordered 
domain state, both in two and three dimensions. This domain state has many of the characteristics of a spin glass, 
as for instance a remanent magnetisation and an irreversibility line scaling like the deAlmeida-Thouless line For 
large disorder the domains have been shown to have a fractal structure []5j with a broad distribution of domain sizes 
and with scaling laws strongly deviating from the usual Imry-Ma assumptions ^ which are thought to be correct in 
the limit of small disorder. 

We focus on the dynamics of the 2D DAFF during relaxation from a long-range ordered system to the fractal-domain 
state and describe this process in terms of a certain growth process dominated by thermal activation. In order to 
prove this approach in detail we use two numerical techniques to investigate 2D systems with a size of 700 x 700 on 
a square lattice and a dilution of 30%. 

In Monte Carlo (MC) simulations we use the standard heat-bath algorithm. Due to the slow dynamics of 
the DAFF at low temperatures it is extremely difficult to investigate equilibrium properties using MC techniques. 
Therefore, additionally we calculate the highly non-trivial ground states of the DAFF exactly using methods known 
from graph theory m . The Ising system is mapped on a graph and the maximum flow through the graph is determined 
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by the Ford-Fulkerson algorithm j9|. The combination of these two methods, MC simulation and exact ground state 
(EGS) calculation, allows a precise investigation of both dynamics and equilibrium properties at finite temperatures 
and at T = 0. 



II. GROUND STATE PROPERTIES 



The geometrical properties of domains in random magnets are important as a starting point for theoretical consid- 
erations as well as for the interpretation of experimental results. In order to investigate the structure of the domain 
state of a simulated DAFF we perform a cluster analysis with a suitable adjusted Hoshcn-Kopclman type algorithm 
p0| . This algorithm pieces the system into domains. A domain is uniquely defined as a group of spins which are 
connected and antiferromagnetically ordered. In this way, it is possible to compute directly the volume v (number 
of spins), surface F (number of unsatisfied bonds) and radius R (root of the mean squared distance of spins) of 
the domains formed (note that in two dimensions v is the area and F the border of the domains). Also one can 
compute the corresponding energy-relevant quantities, like domain wall energy E w (number of broken bonds) and 
volume-magnetization M v for all domains. 

The fractal properties of the domain state calculated by MC simulation have been reported earlier |^,|| . For high 
disorder we discovered a fractal and interpenetrating structure of the domains fulfilling scaling relations which are 
strongly deviating from the assumptions of the Imry-Ma argument ||. These earlier findings are also correct for 
the EGS reported here for the first time for the 2D DAFF. Fig. |l| shows the volume-magnetisation M v (v), i. e. the 
mean volume- magnetisation M v of domains with size v. The scaling behavior is identical with that for the domains 
calculated by MC simulation - which are inevitably out of equilibrium. For 3D systems this has been shown earlier 
||. As Fig. [l] shows the scaling relation M v (v) ~ v dm holds for large enough v. All scaling relations of the DAFF 
obtained so far are summarized in Table Q. 

An important quantity is the distribution of domain sizes found in the domain state of the 2D DAFF. In the next 
section we will use this distribution function for a description of the dynamics of the DAFF. Fig. ^ shows the number 
of domains N for a given volume v from EGS calculations. Since the number of domains is a strongly fluctuating 
quantity the data are averaged over intervals of volumes Aw. The width of these intervals is increased exponentially 
so that we have a constant distance between points on a logarithmic scale. The data are well described by a power-law 
with an exponential cut-off, 

N(v) =N v- dn e- v ' Vo . (2) 

The lines in Fig. |2|are fits to Eq. |[ The cut-off parameter vq is field dependent. Although the fractal dimensions of 
the domains of the DAFF are the same as those known from percolation |ll[] the exponent d n rs 1.5 is between that 
for critical percolation clusters (d n — 187/91) and that for lattice animals (d n = 1). A possible small field dependence 
of this exponent cannot be ruled out. Also, the isolated clusters from the percolation problem lead to the systematic 
deviation from the power law in Fig. ^| for very small volumes v < 5. Since there is no infinite domain in the system 
we do not consider a streched exponential cut-off in Eq. || (the latter is known from the percolation problem if there is 
an infinite cluster). Summed up, the domains of the DAFF have the structure of the lattice animals of the percolation 
problem but with a different exponent d n in the distribution function. 

The distribution function Eq. ^ is a central quantity of a domain state since other quantities follow from this 
distribution. E. g. the order parameter M s per spin can be rewritten as the sum over all domain sizes v of the 
corresponding order parameters times the number of domains of that size, 

oo 

M s = J dv n(v)M s {v), (3) 
1 

where n(v) is the distribution of domain sizes normalized in such a way that M s is the staggered magnetisation per 
spin. 

In the next section we will use the equations above as the starting point for a description of the dynamics of the 
DAFF. 



III. DYNAMICS 



We focus on the dynamics of the two-dimensional DAFF. As was discussed in the preceeding section the equilibrium 
state of the system in a finite field consists of antiferromagnetically ordered domains with a size distribution following 
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Eq. H When we start the dynamics with an antiferromagnetically long-range ordered system how does the system 
evolve into this fractal-domain state? 

Fig. |^ shows a ground state configuration from EGS calculations (bottom) and two configurations of the same 
system (i. e. with same configuration of vacancies) during the MC simulation. For these pictures we restrict ourselves 
to a small system of size 100 x 100. We start the simulation with a long-range ordered state (white) and switch on a 
field B = 1.5 and a finite temperature T — 0.4. 

The ground state of our system is essentially unique (see Ref. || for details). It is a domain pattern that is stored 
in the system through the vacancy configuration and the field and it is the equilibrium state of the system. How is 
this domain pattern reconstructed in our Monte Carlo simulation? 

As Fig. [| suggests, during the relaxation from the long-range ordered initial state to the fractal-domain state the 
dynamics of the system can be understood as a kind of growth process. For shorter times only small domains arise. 
Interestingly, these small domains are located at those places where later the larger domains will be. At later times, 
also larger domains occur and the size of the largest domains increases in time. The pictures lead us to the following 
description of the growth process: the equilibrium distribution of domains develops starting with the smaller domains 
and increasing the maximum domain size in time. 

In order to investigate this process quantitatively we simulate the size distribution of domains during the relaxation 
process. Fig. ^ shows the equilibrium distribution from EGS calculations and distributions of domains that are reversed 
after five different times of MC simulation. For a given time t the system has obviously a distribution that is equal 
to the equilibrium distribution n(v) corresponding to Eq.|^ up to the largest domain size v m (t) that has developed 
at that time by thermal activation. Note that this time dependent cut-off in the distribution function is very abrupt 
since Fig. ||] shows data over a range of eight decades. These considerations lead to the following description of the 
time development of the order parameter M s . We start with a long-range ordered system, i. e. M s (t — 0) = M si w 1. 
[j] After switching on the field the staggered magnetisation decreases to M s (t — > oo) = M soo w (except of finite-size 
effects) due to the growth of domains of the other phase: 

v m (t) 

M s (t) = M si +2 J dvn(v)M s (v) (4) 
l 

Since a domain is defined as a group of spins which are antiferromagnetically ordered it is per definition M s (v) = —v 
for the reversed domains so that 



v m (t) 



M s (t) = M sl -2 J dv n v 



1— dng—v/vo 



M soo + 2n / dv y^^e-v/vo 



v m (t) 

= M soo +j{v m (t)) (5) 
where the last line defines the function "f{v m (t))- no is defined by the normalisation condition 

2 / dv n v 1 - d "e- v / V0 + M soc = M m . (6) 



According to the preceeding section the exponent is d n — 1.5 and the field dependent cut-off parameter can be obtained 
from EGS calculations. Thus j(v m (t)) is well defined without any unknown parameters. Since the integration cannot 
be done exactly it will be calculated numerically. Furthermore, the integral j(v m (t)) is related to the so-called 
incomplete Gamma- function 



1 Note that linit^o M s < 1 since due to the dilution there are spins or clusters of spins that are not or only weakly connected 
to the infinite cluster. These clusters may have finite magnetisations which follow the external field immediately since there 
are no relevant energy barriers. 
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T(a,x) = / dyy^e-y (7) 



the properties of which we will discuss later. 

It is highly plausible that the largest existing domain size v m (t) for a time t is connected to a thermal activation 
energy AE(t) = Tln(t/r) according to 

AE(t) = E b vt(t) (8) 

where db is an exponent connecting the size of domains with a free-energy barrier AE which has to be overcome 
during the build-up of a domain of size v m . This assumption of thermal activation gives the time dependence of M s 
in Eq.| 

In order to prove the validity of our approach we use MC simulation data for the decay of the order parameter. 
In Fig. |^ we show data for three different values of the magnetic field B — 1.5,2.0,2.5. The way we evaluate the 
simulation data is the following: 

i) First of all, due to the equations above data for the relaxation for different temperatures should collapse in a 
scaling plot M s versus any function /(Tln(t/r)) where r should be a microscopic time scale for the growth of the 
smallest domains. In Fig. || the data for four (B — 1.5) respectively three (B = 2.0,2.5) different temperatures and 
a time window ranging from 95 MCS to 200000 MCS are seen to collapse very well using r = 3 MCS. We do not 
consider time scales much shorter than 100 MCS because then thermal fluctuations cannot be neglected. 

ii) Finding the correct barrier exponent d b and the prefactor Eb we can determine v m (t) = ((T/E b ) ln(i/r)) 1/db and 
following Eq. || the data for M s (t) — M soo should be described by j(v m (t)). As Fig. || shows this works perfectly for 
the complete range of the scaling variable v and for three different fields using d b = 0.25. Note that the reason for 
the field dependence of -f(v m (t)) is the field dependence of the cut-off parameter u which has been obtained by EGS 
calculations (see Fig. ^]). 

The data for the two higher fields are less accurate since in this case the dynamics is faster and has a shorter 
observation time. Additionally for the higher fields we restricted ourselves to smaller lattices (400 x 400 instead of 
700 x 700 for B = 1.5) since the EGS calculations are very time consuming. However, no field dependence of the barrier 
exponent can be seen. The prefactor E b is slightly field dependent (E b = 0.60(5 = 1.5), 0.74(£ = 2.0), 0.67(5 = 2.5)). 
Note that since the barrier exponent is very small we are close to the case AE{t) ~ \n(v m (t)) which is the limit d b — > 
in Eq. § 

Of course, other quantities can be described in a similar manner. E. g. for large enough volumes v (see Fig. |l|) the 
magnetisation M can be written correspondingly as 

M(t) = Mi + 2 [ dv n M v d ™- d »e- v / Vo (9) 



The magnetisation M is a quantity that can be measured experimentally. Its precise behavior has to be obtained 
from a numerical integration of Eq. [)| Two for the description of experimental data relevant limits can however be 
obtained asymptotically as limits of the Gamma-function. 

The first limit is that of short times and small temperatures where the exponential cut-off is irrelevant which means 
((T/E b )ln(t/T)) 1 ' db = v <C i>o and consequently e~ v / v ° « 1 resulting in 

M(t) w Mi + 2n M (v(t)) dm - d " +1 

= Mi + 2noM Q (— ln(t/r)) d '""4" +1 . (10) 

This is a logarithmic law with a temperature-independent exponent obtained from the fractality and the distribution 
of the equilibrium domains and from the free-energy barriers that have to be overcome to build up these domains. 
Similar logarithmic laws are usually used to describe experimental data of the dynamics of the DAFF (for a review 
of the experimental work see fl2)). Note that Mi is not zero but has a small finite value for the same reason that has 
been explained earlier - see the discussion of Af S j. 

The second limit is that close to equilibrium, i. e. for long times and higher temperatures. Here we can use an 
asymptotic form of the Gamma-function 



lim T(a,x) ~ x a e 
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The restriction to the most relevant term - the e-function - results in 



^ t -T/(E bVo ) for dt = 1 



(11) 

(12) 



The generalized power law Eq. |TT] has been used to fit experimental data as well as data from MC simulations for the 
remanent magnetisation of the 3D DAFF For db = 1 it simplifies to a simple power law with the exponent 

proportional to the temperature. This law has also been discussed as a description of the remanent magnetisation of 
the 3D DAFF and it is typical for the decay of the remanent magnetisation of spin glasses (see |L6| for a review). 



IV. CONCLUSIONS 



We derived a simple description for the dynamics of the 2D DAFF during the relaxation from long-range order to a 
fractal-domain state. The basic fundaments of this description are the fractality and the distribution of the domains 
the domain state of the system consists of in the limit of high disorder. This is the main extension to earlier theories 
for the case of weak disorder based on scaling assumptions for domains which are compact and non-fractal with for a 
certain time scale unique relevant length scale - the domain radius (see e. g. [|||l7j for scaling theories of random-field 
models and fll7| , p"8| for a description of domain dynamics in terms of domain wall motion). 

As a result of our approach which is based on measured quantities of the domain state the dynamics can be 
described with the assumption of thermal activation analytically and is given by an incomplete Gamma-function. As 
one advantage this function connects naturally two for disordered systems often observed laws as limiting cases — 
the logarithmic law in the limit of short times and low temperatures and the (generalized) power law in the limit of 
late times and higher temperatures. Apart from that our approach illustrates that from the fact that a system has 
a distribution of domain sizes instead of one dominating size follows that the mean domain size of that distribution 
can grow following a power law (Eq. |l2] ) even if the largest domain in the system is build by thermal activation, i. e. 
within a logarithmic time scale (see e. g. fll6[ for a summary of this apparent contradiction in spin glasses). 

The agreement with simulations is excellent. However, so far we are not aware of any direct experimental observation 
of this kind of dynamics in 2D DAFF although as mentioned above similar laws are often observed in the dynamics of 
disordered magnets and we assume that a similar description might also work for other strongly disordered systems 
like random-field systems or spin glasses. 
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FIG. 1. Magnetisation M v versus volume v of the domains of the fractal- domain state of the 2D DAFF. Comparison of MC 
simulations (avarage over 50 systems of size 199 x 198) and ECS calculations (avarage over 10 systems of size 400 x 400). The 
latter curve is shifted by a factor 2. B = 1.5. All quantities are dimensionless. 

FIG. 2. N versus v from EGS calculations for different fields B. Apart from the average over 10 systems of size 400 x 400 the 
data are additionally averaged in such a way that they show the relative number of domains having a volume within a certain 
intervall Av. The lines are fits to Eq. ^. 

FIG. 3. EGS of a 100 x 100 system (bottom) and two configurations of the system during the MC simulation of the relaxation 
process after roughly 500 MCS (top) and 130000 MCS (center). T = 0.4, B — 1.5. The two antiferromagnetic phases are 
represented in black and white, the vacancies are grey. 

FIG. 4. N versus v during MC simulation of the relaxation process and for the EGS of the same 700 x 700 system (see Fig. ^ 
for explanation). T = 0.4, B = 1.5. 

FIG. 5. Scaling plot for the order parameter AM, = M a (t) — M soo . Data are shown for three different fields, B = 1.5, 2.0, 2.5 
(from above) and four (B = 1.5) respectively three (B = 2.0,2.5) different temperatures. The time ranges from 95 MCS to 
200000 MCS. The solid lines are the theoretical curves following Eq. | with v = 9000(B = 1.5), 1000(B = 2.0),400(B = 2.5). 



TABLE I. Scaling relations for the DAFF 
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